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Abstract — We  consider  a  formation  control  problem 
for  a  robotic  network  with  limited  communication  and 
controlled  motion  abilities.  We  propose  a  novel  control 
structure  that  organizes  the  robots  in  concentric  layers 
and  that  associates  to  each  layer  a  local  leader.  Through  a 
load  balancing  algorithm  on  the  asynchronous  network 
of  layers  we  allocate  the  desired  number  of  robots  on 
each  layer.  A  final  uniform  spreading  algorithm  leads  the 
robots  to  a  lattice-like  formation.  This  novel  distributed 
communication  and  control  algorithm  runs  in  linear  time 
in  the  worst  case. 

I.  Introduction 

Models  for  robotic  networks  are  the  subject  of  the 
early  work  in  [1],  where  some  impossibility  results 
are  established  for  certain  formation  control  problems. 
More  recently,  Martinez  et  al  in  [2]  propose  a  detailed 
model  and  analyze  the  time  complexity  of  basic  ren¬ 
dezvous  and  deployment  algorithms.  For  many  of  the 
resulting  linear  dynamical  systems,  the  worst  case  time 
complexity  is  of  order  @(n2logn)  (rendezvous)  and 
()(n:’>  log  n)  (deployment).  Sharma  et  al  [3]  characterize 
the  time  complexity  of  sensor-based  vehicle  routing 
in  2  dimensions,  i.e.,  of  a  distributed  point  to  point 
motion  planning  problem.  In  their  formulation,  size  and 
speed  of  the  robots  decrease  with  an  increasing  size  of 
the  network;  in  summary  the  worst  case  complexity  of 
the  vehicle  routing  problem  is  0(n)  and  the  expected 
time  complexity  if  Q(^/n).  Smith  and  Bullo  in  [4] 
consider  a  target  assignment  problem  and  characterize 
its  time  complexity  as  being  @(n)  in  the  worst  case.  It 
is  not  completely  fair  to  compare  all  these  results  as  the 
models  they  adopt  have  small  but  important  differences. 
Additional  related  work  includes  the  study  of  scaling 
laws  in  statistical  mechanics,  percolation  theory,  and 
wireless  communication;  e.g.,  see  [5],  [6],  [7],  [8]. 

The  contributions  of  this  paper  are  twofold.  We 
consider  a  formation  control  problem  for  a  robotic 
network  with  appropriate  scaling  laws  inspired  by  the 
thermodynamic  limit  in  statistical  mechanics  and  by 


the  memory  and  communication  rate  restrictions  typ¬ 
ical  in  distributed  algorithms.  Our  formation  control 
problem  amounts  to  the  problem  of  steering  robots  to 
“contiguous”  vertices  of  a  lattice  formation  in  the  plane. 
Our  first  contribution  is  to  show  that  this  problem  has 
worst  case  time  complexity  of  order  0(n ).  The  second 
contribution  is  the  explicit  design  of  an  communication 
and  control  law  that  achieves  this  upper  bound.  Our 
novel  approach  is  based  on  the  distributed  computation 
and  representation  of  a  control  structure  called  “rotating 
layers.”  The  key  idea  is  to  organize  the  robots  in 
concentric  layers  and  to  orchestrate  their  motion  in 
such  a  way  that  the  resulting  network  is  connected  over 
certain  intervals.  The  main  algorithmic  contribution  is 
a  load  balancing  algorithm  for  an  asynchronous  chain 
network  with  bounded  time  complexity. 

The  paper  is  organized  as  follows.  Section  II  presents 
the  formal  model  with  a  particular  attention  to  appro¬ 
priate  scaling  laws  and  complexity  notions.  Section  III 
introduces  our  novel  control  structure  and  establishes 
its  properties.  Section  IV  contain  the  main  algorithmic 
contribution,  namely  a  simple  load  balancing  algorithm 
and  its  analysis.  Finally,  Section  V  describes  simple 
strategies  on  how  to  place  the  robots  in  their  desired 
positions  inside  the  control  structure. 

II.  Modeling  the  robotic  network 

Let  n  be  the  number  of  robots  in  the  network,  each 
with  a  unique  identifier  (UID).  Let  /  =  {1, . . .  ,  n}  be 
the  set  of  unique  identifiers.  We  assume  continuous¬ 
time  motion  and  communication  at  discrete  times. 
Each  robot  has  the  following  sensing,  computation, 
communication,  and  motion  control  capabilities.  The 
memory  and  transmission  capability  of  the  processor 
is  0(log(n))  bits.  The  zth  robot  occupies  a  location  p* 
in  a  compact  region  Qcl2,  and  it  moves  according 
to  the  continuous-time  control  system 

Pi  =  Ui, 


with  | it |  <  "Umax-  We  shall  let  the  robots  to  always  move 
at  the  maximum  speed  ttmax-  The  sensing  and  com¬ 
munication  model  is  the  following.  The  processor  of 
each  robot  can  sense  its  position  with  1  Jyfn  accuracy. 
Consequently,  each  robot  can  represent  its  position  with 
0(log(n))  bits.  Every  robot  is  capable  of  transmitting 
its  position  to  any  other  robot  within  a  closed  disk  of 
radius  rcomm  £  R+.  We  assume  that  rcomm  and  umax 
scales  as  1/ y/n  because  of  congestion.  The  inspiration 
for  this  scaling  comes  from  the  thermodynamic  limit 
in  statistical  mechanics. 

The  n  robots  are  initially  randomly  placed  in  Q. 
The  communication  graph  is  the  rcomm-disk  graph, 
i.e.,  the  graph  Cdisk  with  vertices  pi,...,pn  and  edges 
( Pi,Pj )  for  all  i  /  j  if  and  only  if  \\pt  —  Pj\\  < 
rcomm.  We  assume  that  the  initial  placement  of  the 
robots  pi(0), . . .  ,pn( 0)  is  such  that  their  communica¬ 
tion  graph  C/disk  (pi(0), ...  ,pn(0)  is  connected. 

The  coordination  objective  is  as  follows:  place  the 
robots  at  “contiguous”  locations  in  a  hexagonal  lattice 
with  desired  inter-robot  distance  rdes  £  R+  (refer 
Fig  1),  where  rdes  scales  as  1  /y/n  because  of  conges¬ 
tion.  We  assume  that  rc|es  <  rcomm/2. 


Fig.  1.  Hexagonal  lattice  formation. 


III.  A  CONTROL  STRUCTURE  BASED  ON  ROTATING 
LAYERS 

In  this  section  we  outline  the  details  of  a  control 
structure  that  the  robots  compute  as  part  of  their 
initialization  before  executing  the  main  algorithm.  The 
communication  graph  used  for  this  initialization  pro¬ 
cess  is  the  rcomm-disk  graph,  i.e.,  two  robots  are  able  to 
communicate  to  each  other  only  when  they  are  within 
'/"comm  distance  of  each  other.  The  control  structure  is 
constructed  via  a  few  intermediate  steps. 

Step  1:  Central  leader  election 

The  robots  estimate  the  boundaries  of  the  bounding 
box,  i.e.,  the  smallest  rectangle  that  is  aligned  with  the 


Fig.  3.  Rotating  Layers. 


coordinate  axes  and  that  contains  all  the  robots.  The 
robots  then  calculate  their  respective  distances  from 
the  geometric  center  of  this  rectangle  and  the  robot 
(o  that  is  closest  to  the  center  is  elected  as  the  leader. 
A  flooding  protocol  to  compute  the  bounding  box  is 
discussed  in  [9];  it  is  immediate  to  see  that  its  time 
complexity  is  0(n),  where  n  is  the  number  of  robots. 
As  is  well  know  [10],  [11],  the  time  complexity  of 
leader  election  based  on  a  flooding  algorithm  is  also  in 
0(n). 

Step  2:  Rotating  layer  formation 

The  second  step  in  computing  the  control  structure 
is  the  formation  of  rotating  layers  centered  at  the  leader 
Iq.  To  do  this,  the  robots  first  determine  the  layer  to 
which  they  belong  based  on  their  distances  from  the 


leader  1$  (refer  Figure  2).  The  layer  number  of  robot 
i,  denoted  by  L(i),  is  given  by 


m 


~jj  Pi-Plo 

r'dcs 


Because  the  environment  is  bounded,  so  is  the  maxi¬ 
mum  layer  number,  that  we  denote  with  iViaycr.  As  a 
consequence  of  the  connectivity  assumption  and  of  the 
relationship  between  rcomm  and  rc|cs,  one  can  see  that  if 
the  annulus  with  radii  [jrdes,  (j  +  l)^des]  is  empty,  then 
so  are  all  other  annuli  with  larger  radii. 

The  robots  then  move  radially  away  from  pp  to  a 
distance  from  pp  determined  by  their  layer  number. 
The  nominal  distance  from  £q  for  robot  i  is  rnom  (i) 
defined  by 

^nom(i)  —  ^(i)^des' 


Once  the  robots  reach  their  nominal  distances  as  shown 
in  Figure  3,  the  ones  in  the  odd  numbered  layers  rotate 
in  the  clockwise  direction  and  the  ones  in  the  even  num¬ 
bered  layers  rotate  in  the  counterclockwise  direction. 
We  shall  refer  to  such  layers  as  rotating  layers  (refer 
Figure  3).  Without  loss  of  generality,  we  assume  that 
the  speed  umax  is  sufficiently  small,  as  compared  with 
rComm,  so  that  any  robot  in  layer  number  k  is  within 
communication  range  of  all  robots  in  layer  numbers 
k  —  1  and  k  +  1  at  least  once  every  2tt  A'rdcs  /'Umax  time 
intervals.  The  resulting  communication  graph  over  an 
appropriate  time  interval  is  shown  in  Figure  4.  This 
step  takes  0(1)  time. 


Fig.  4.  Resulting  “over-time”  communication  graph 


Step  3:  Layer  leader  elections 

After  all  robots  reach  their  nominal  distances  and 
begin  rotating,  the  following  protocol  is  executed.  The 
objective  is  to  elect  a  leader  l3  for  each  rotating  layer 
j  and  to  count  the  number  of  robots  present  in  each 
layer. 

1:  set  j  :=  1 


2:  robot  £j_\  counts  the  number  of  robots 
in  layer  j,  denoted  by  n\ayer(j),  during 
the  next  27rjrdes/umax  time  interval 
3:  if  niayer(j)  =  0  then 
4:  exit 

5:  else 

6:  among  all  the  robots  present  in  layer 

j,  robot  £j-\  elects  the  one  with  the 
highest  UID  as  lj  and  communicates 

^layer(j)  to  it 

7:  end  if 

8:  Set  j  :=  j  +  1  and  go  back  to  2  : 

For  layer  j,  the  counting  and  leader  election  process 
requires  at  most  47rjVdes/Rmax  time.  The  following 
lemma  characterizes  the  time  complexity  of  the  process 
required  for  the  layer  leader  election  step. 

Lemma  3.1:  Starting  from  the  rotating  layers,  the 
time  required  to  elect  layer  leaders  belongs  to  0(n). 

Proof:  Circumference  of  the  jth  layer  is 

27T7Vdes(ro)  ~  j/y/n-  Robot  speed  is  umax(n)  ~  1  jyfn. 
Flence,  steps  2  :  and  6:  require  O(j).  But  the  number 
of  layers,  Ajaycr  S  Oiffi).  Therefore  the  total  time 
required  to  elect  layer  leaders  is  3  F  O(n).  ■ 

We  are  now  ready  to  state  the  time  complexity 
required  to  form  the  rotating  layer  structure. 

Lemma  3.2:  The  time  required  to  form  the  rotating 
layer  structure  from  an  initially  connected  set  of  n 
robots  belongs  to  0(n). 

Analysis  of  communication  across  and  inside  layers 

As  the  control  structure  is  in  place  and  robots 
rotate  according  to  their  layer  identifier,  the  commu¬ 
nication  graph  across  robots  changes  with  time.  In 
other  words,  if  the  robot  trajectories  are  denoted  by 
(pi(f), . . . ,  pn  (t)),  then  the  time-dependent  communi¬ 
cation  graph  is  t  i->  Gdisk(pi(t), . . .  ,pn(t)).  Because 
the  number  of  layers  belongs  to  0(^/n),  this  time- 
dependent  graph  has  the  following  useful  properties 
(refer  Figure  5): 

(i)  for  all  j  >  0,  there  is  an  edge  between  l3  and 
any  robot  in  layer  j  +  1  at  least  once  every 
2vrjrdes/Rmax  G  0(y/n)  time  instants,  and 

(ii)  each  leader  £j  can  broadcast  a  message  to  all 
robots  in  its  own  layer  j  in  time  47rjrdes/Rmax  E 
0(y/n). 


IV.  Layer  balancing 

The  hexagonal  lattice  can  be  viewed  as  a  layer 
formation  around  a  central  point  with  6 j  points  in  layer 
number  j.  In  order  to  arrange  the  robots  into  such 


Transmission  delay  0(j) 


Fig.  5.  Asynchronous  chain  network 

a  configuration,  we  first  need  to  distribute  the  robots 
among  the  layers  so  that  we  have  6 j  robots  in  the 
jth  layer.  We  want  this  property  to  hold  for  as  many 
consecutive  layers  as  possible  starting  with  layer  1.  We 
call  this  the  layer  balancing  problem.  To  solve  this 
problem,  we  introduce  a  useful  finite  state  machine  £ 
and  design  a  load  balancing  algorithm  for  £.  We  later 
show  how  to  implement  this  algorithm  on  the  robotic 
network. 

The  state  of  the  finite  state  machine  £  is  the  number 
of  robots  in  each  layer.  Let  N*  be  the  smallest  positive 
integer  such  that  N*(N*  +  1)  >  n/3  (this  will  be  the 
final  number  of  balanced  layers).  The  number  of  states 
of  £,  denoted  by  ATdyn,  is  max{A}ayer,  A'*}.  Hence, 
the  state  of  £  can  be  described  by  the  tuple  x{t)  = 
(xi(t), . . .  ,xNiiyn(t)),  where  xj( 0)  =  n\ayer(j)  VI  < 

j  <  Allayer  and  Xj( 0)  =  0  V Allayer  +  1  <  j  < 

Ndyn.  In  addition  to  this,  we  define  auxiliary  state 
variables  y(t)  =  (yi(t), . . . ,  yjvdyn(*))  and  z(t)  = 
(z\(t), . . . ,  ZNdyn(t)).  These  variables  are  initialized  as 
follows:  yi(0)  =  1,  yj{ 0)  =  0  Mi  >  1  and  zy (0)  = 
0  Mi  >  1.  These  auxiliary  variables  y(t)  and  z(t)  take 
values  in  {0, 1}  and  their  values  influence  the  evolution 
of  £.  The  evolution  of  £  obeys  the  following  discrete¬ 
time  linear  system: 


'+1' 

'  0  ' 

-1 

0 

Ul  +  ... 

0 

+1 

0 

-1 

where  Uj  £  Z  is  a  control  signal  establishing  how 
many  robots  should  transfer  from  layer  j  +  1  to  j. 
For  the  system  £,  the  layer  balancing  problem  amounts 
to  the  selection  of  an  input  trajectory  that  steers  x(t) 
from  the  initial  state  (xi(0), . . .  ,xjvdy„(0))  to  a  final 
state  (xi(T), . . . , xNiyn(T)),  where  Xj(T)  =  6 j  VI  < 

j  <  N*  -  1,  xN*  =  n  -  J2jLixj(T)  and  xj(T )  = 
0  MN*  +  1  <  j  <  Ndyn  for  some  T  >  0.  This  is  akin 
to  a  load  balancing  problem  on  an  asynchronous  chain 
network. 

We  assume  that  only  one  control  signal  can  be 


nonvanishing  at  each  instant  of  time.  In  the  following 
algorithm  we  denote  by  £  <—  (a,  (3)  the  transition  of  £ 
generated  by  the  control  signal  ua  =  (3  with  all  other 
inputs  equal  to  zero. 

Incremental  balancing  algorithm 
1:  for  j  :=  1  to  (ATdyn  -  1)  do 
2:  if  (yj  =  1)  AND  (xj  >  6 j)  then 

3:  {push  excess  robots  to  next  layer} 

4:  E<-  (j,6j  -Xj) 

5:  {relinquish  the  token  to  next  leader} 

6:  set  yj+ 1  :=  1,  yj  :=  0  and  Zj  :=  1 

7:  end  if 

8:  if  (yj  =  1)  AND  (xj  >  6 j)  then 

9:  {push  excess  robots  to  previous 

layer} 

10:  £  <-  (j  -  l,Xj  -  6j) 

11:  end  if 

Lemma  4.1:  The  incremental  balancing  algorithm 
solves  the  layer  balancing  problem  in  time  0(n)  in 
the  worst  case. 

Proof:  We  only  provide  a  sketch  of  the  proof  here 
and  refer  the  reader  to  a  forthcoming  manuscript.  If 
no  layer  leader  ever  waits,  then  leader  j  requires  time 
of  order  O(j)  to  execute  the  algorithm  and  the  time 
complexity  is  3  =  0(n).  On  the  other  hand, 

if  the  layer  leader  i  waits  for  robots  to  arrive  from  the 
layer  j  with  j  >  i,  then  one  can  deduce  that  the  waiting 
time  is  at  most  0(J2l=i  k)  and  that,  after  i  sets  its  Z{ 
variable  to  1,  then  no  more  layer  leader  will  wait  until 
layer  leader  j  at  least.  In  summary,  the  total  waiting 
time  is  at  most  j)  =  0(n).  ■ 

Next,  we  explain  how  the  robotic  network  equipped 
with  the  rotating  layer  structure  can  implement  the  IN¬ 
CREMENTAL  balancing  ALGORITHM  in  a  sequential 
asynchronous  manner. 

We  begin  by  designing  some  communication  mes¬ 
sages  between  layer  leaders  that  encode  various  inter¬ 
layer  cooperation  requests.  Based  on  the  following 
requests  one  can  easily  envision  how  a  distributed 
algorithm  for  the  synchronous  robotic  network  can 
be  designed  to  exactly  implement  the  INCREMEN¬ 
TAL  BALANCING  ALGORITHM  for  the  asynchronous 
chained  network. 

Relinquish  from  j  to  j  +  1:  Provided  layer  j  +  1 
exists,  leader  ij  communicates  to  Zj+ 1  that  it  is  now 
its  turn  to  evaluate  if  the  layer  j  +  1  is  balanced  or  not. 

Push  excess  robots  from  j  to  j  + 1:  Provided  layer 
j  contains  at  least  6 j  robots,  leader  ij  requests  the 
excess  robots  in  its  layer  to  move  to  layer  j+1  (possibly 


creating  it  and  subsequently  electing  a  leader).  The 
counters  on  layers  j  and  j  +  1  are  updated  through 
a  communication  between  tj  and  ^j+ 1- 

Push  excess  robots  from  j  to  j  —  1:  Provided  layer 
j  contains  at  least  6 j  robots,  leader  £3  requests  the 
excess  robots  in  its  layer  to  move  to  layer  j  —  1.  The 
counters  on  layers  j  and  j  —  1  are  updated  through  a 
communication  between  tj  and  £j-\. 

V.  Uniform  spreading 

This  phase  is  executed  sequentially  starting  from 
innermost  layer  and  progressing  outwards.  We  provide 
an  informal  description  in  three  steps: 

(i)  When  the  leader  t3  is  in  charge  of  the  uniform 
spreading  of  the  layer  j  + 1,  starting  from  £j+\,  it 
does  the  following  two  things  as  it  encounters  a 
robot  i  in  layer  j+ 1.  First,  it  determines  the  order 
number  of  the  robot  i  in  the  layer,  denoted  by 
order (i),  i.e.,  1  for  the  leader  ^j+ 1>  ^  for  the  next 
robot  and  so  on.  Second,  it  transmits  the  STOP 
signal  and  the  quantity  9(i)  =  27r(i— l)/niayer(j+ 
1)  to  the  robot  i. 

(ii)  When  a  robot  i  in  layer  j  +  1  receives  these 
messages  from  the  leader  tj,  it  makes  one  full 
rotation  and  then  stops  at  a  position  which  ex¬ 
tends  an  angle  6(i)  at  £q  with  respect  to  the 
north  direction.  (One  way  to  do  this  is  storing 
coordinates  of  one  point  which  is  exactly  north 
of  4). 

(iii)  After  stopping,  t3+\  takes  charge  of  the  uniform 
distribution  of  the  layer  j  +  2.  All  the  other 
robots  in  the  layer  j  +  1,  move  inwards  by  an 
appropriate  distance  until  they  are  at  a  vertex  on 
the  hexagonal  lattice.  (One  way  to  do  this  is  by 
making  every  robot  move  radially  inwards  until 
it  cannot  move  further  without  being  closer  than 
rdes  to  any  other  robot). 

On  the  completion  of  this  uniform  spreading  step, 
all  the  robots  have  occupied  the  vertices  of  a  regular 
lattice.  This  step  takes  0(n)  time. 

VI.  Summary 

In  summary,  we  can  now  state  a  result  for  the  total 
time  required  to  solve  the  lattice  formation  problem 
using  the  procedure  in  this  paper. 

Theorem  6.1:  The  time  required  to  solve  the  lattice 
formation  problem  for  a  robotic  network  described  in 
Section  II  belongs  to  0{n).  This  bound  is  achieved 
by  the  sequential  computation  of  the  rotating  layers 


and  execution  of  the  incremental  balancing  and  uniform 
spreading  algorithms. 

In  future  work  we  envision  studying  worst-case 
lower  bounds  as  well  as  investigating  stochastic  sce¬ 
narios  in  which  the  initial  robot  positions  are  randomly 
uniformly  placed. 
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